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Abstract 

We study non-superposition effects in the Dirichlet Casimir interaction 
energy for N boundaries in d spatial dimensions, quantifying its de- 
parture from the case of an interaction where a superposition principle 
is valid. We first derive some general results about those effects, and 
then show that they only become negligible when the distances be- 
tween surfaces are larger than the sizes of each individual surface. We 
consider different examples in one, two and three spatial dimensions. 

1 Introduction 

Among the many interesting and distinctive features of the Casimir effect [T], 
not the least important is the fact that the forces it produces between (more 
than two) conducting surfaces do not satisfy, in general, a superposition 
principle. In other words, when dealing with more than two objects, the 
interaction energy cannot be written as the sum of the interaction energies 
corresponding to all the possible object pairs. As a consequence, knowledge 
of the energy of a system before the addition of an extra surface may seem to 
be of little help, since there is no obvious way to include the extra interaction 
terms. Not unexpectedly, this non-superposition property is shared by the 
van der Waals interaction [2]. 

Quite apart from its theoretical interest, it should be noted that this kind 
of phenomenon may also be of practical relevance, since it could be helpful in 
some approximate calculation schemes, in situations where nonlinear effects 
are small. For example, if there exists a regime where superposition is approx- 
imately valid, one should expect the dominant term in the Casimir energy to 
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be akin to a two body interaction potential, albeit with a non-Coulombian 
potential, plus corrections. Under some assumptions, those corrections can 
be, as we shall see, small perturbations. 

In this paper, we first investigate the non-superposition effects in a quite 
general approach. We then argue that when those effects are small, a pertur- 
bative expansion naturally suggests itself. We then discuss and apply that 
approximation within the context of different examples. 

This paper is organized as follows: in section [2] we first review the func- 
tional approach to the calculation of the Casimir energy. In [3l we deal with 
the study of non-superposition effects, relating them to a perturbative ex- 
pansion inSlwith examples in 1 + 1, 2 + 1, and 3 + 1 dimensions. 
In [HI we present our conclusions. 

2 The method 

In order to analyze properties of the Casimir energy, it is convenient to 
introduce one of its concrete representations. We shall use here one that is 
based on the functional integral formalism introduced in [3l[4l[5]. 
What follows is a review of its main aspects, adapted to the problem at hand. 

Denoting by the Euclidean vacuum amplitude for a real mass- 

less scalar field Lp in the presence of Dirichlet surfaces E^*^), a = 1, . . . , A^, 
the total vacuum energy Eq may be written as follows: 

where T is the extent of the (imaginary) time interval, and Zq is the free (no 
surfaces) vacuum amplitude. The role of the latter is just to fix the vacuum 
energy to zero when there are no surfaces. 

On the other hand, the vacuum amplitude may be written as a functional 
integrafl: 

ZfjS^'^)}] = j [V^] e-^°(^) , (2) 

where 5*0 is the free Euclidean action, which for a massless scalar field in 
d + 1 dimensions reads: Sq{(p) = | / dfi^^x {dipY, and [Dip] denotes the 
path integral measure corresponding to a scalar field which satisfies Dirichlet 
boundary conditions on each surface S^"^ 

^In our use of the functional integral formalism, we follow the approach and conventions 
of [6]. 
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It is quite useful to write that measure in the equivalent way: 

N 

[V^] = X II [^] , (3) 

a=l 

where we introduced a (5-functional for the field on each surface. 

In what follows, we focus on the d = 3 case, although every step will have 
its analogue for different numbers of dimensions. The changes to the final 
expressions required to deal with d 3 are described in [3l 

Thus, assuming that {a^ja"^) y^"'\o-) (y*-"^ G R.^^^) is a parametriza- 
tion of Tj^"'\ we introduce an auxiliary field (^'•"^(r, cr) to exponentiate each 
functional delta: 

(no sum over a), where g^°'\a) = det [^'^'^^(cr)] , (a,/? = 1,2) is the determi- 
nant of the induced metric g^^^ on the surface, and r = xq. In terms of the 
previous parametrization, 

9 = —Q^a (no sum over a) . (o) 

Equation Q guarantees the (necessary) reparametrization invariance on each 
surface, assuming that the auxiliary fields behave as scalars under those 
transformations. 

Inserting ^ into ([3]), we are left with an equivalent functional integral 
expression for Z[{S*^"^}]: 

= j (^n^?^"^) / exp{-5o(¥^) + I I d'xJ{x)^{x)] 

(6) 

where we introduced J{x) = J2a=i J^"'K^)i with: 

J('^) (x) = J drd^a ^Jg^-\a) ^^^^ (r, a) 6{xq - r) 5^^^ [x - y ("^^ (a)] . (7) 

Performing now the (Gaussian) integral over the field, the result may 
be put, in a condensed form, as follows: 

= Zo X [v^ e-^^«) (8) 
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where = and = / e-^(^«). 

In ([HI), T)^ denotes the integration measure for all the auxiliary fields (we 
assume there is more than one boundary) and S's is a 'nonlocal action' for 
those fields: 

, „ „ TV 

X Mi^a,)(T,o-T\a')e\T\a')\ , (9) 

where each matrix elements of Ai may be expressed in terms of the scalar 
field propagator /C: 

M^ab){r,a;T',a') = ^g(-){a) /C(t - r'; y^^\a) - y<^'\a')) g^'){a') , (10) 
which, in (i = 3, may be written as follows: 



/C(xo,x) = I ^e^-^«/C(c.,x) 



Taking advantage of the time independence of the physical system considered, 
we Fourier transform in time the auxiliary fields, to write their action in the 
following way: 

AT 



J J J a,b=l 

X M^ab){u;;a,a')i'-'\u,a')] , (12) 

where the tilde on the fields denotes their corresponding Fourier transformed 
versions, and: 

Miab){uJ;(J,a') = ^gi-){a)IC^ab){uJ;a,a')^g(b){a') , (13) 

where 

/C(,,)(a;;a,a') = IC[uj;y^''\a) - y^'\a')] . (14) 
Since the integral over the auxiliary fields is Gaussian, we have: 

^ = {det[M(^ab){u;-a,a')6{u-u;')]] \ (15) 



where the determinant refers to the continuous indices cr, cr', as well as 

the discrete ones a, b. 

Then, recalling the relation between the vacuum functional and the vac- 
uum energy Eq we find for the latter the expression: 



where the trace affects the a, a' and a, h indices (the trace over frequencies 
has been explicitly dealt with by means of the integral). We shall reserve the 
symbol 'tr' for the cases where a trace over just the (continuous) cr, a' indices 
is needed. 

Note that no subtraction of the would-be Casimir 'self-energies' has yet 
been performed; this step will be considered in the next section. 

3 Non superposition 

Since we are interested in the Casimir interaction energy, we will first ex- 
tract the self-energies of the surfaces. Besides, those energies are additive 
quantities, insensitive to the phenomenon we wish to consider. 

That extraction can be done by factorizing a diagonal matrix constructed 
from the a = h elements of A^: 



Mi^ab){oo] a, a') = / dV" ^ P(ac)(^; ^, o") T^,t)iuj; a", a') (17) 




(16) 



N 



c=l 



where 



(no sum over a). By construction, Tioj) has the matrix elements: 



(18) 





(19) 



(no sum over a). This factorization implies that det(7Vl) = det(P) det(T); 
thus, recalling (fTSl) and (fT6l) . we may write: 



N 




(20) 



a=l 



where 




(21) 
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is the Casimir self-energy of the object labelled by the index a, and: 











/ ZTT 





[22] 



As already advanced, the self-energies, as seen from ([20l) . are additive. Be- 
sides, they do not contribute to the Casimir forces between the surfaces 0, 
since they are independent of their relative distances. ^ ^ 

Let us then consider the interaction term, Ei. It depends on T = X + T' 
where X is the identity matrix (in both discrete and continuous indices) and 
T' has vanishing diagonal (a = h) elements. Moreover, for a 7^ 6, it coincides 

with f(,fe) of (na). 

As a final step to obtain our main result, we derive a different (but equiv- 
alent) expression for 7(^6), such that formula for the interaction energy does 
not contain explicit factors of the metric. 

To that end, we introduce G^'^\uj] a, a'), the inverse of /C(aa)(a;; a, a'): 



Then: 



and: 



a — a 
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where 
0(afe)(a;; a, a') 



/rfV"G('^)(a;;a,(T")/C("'')(cu;a",a') if a 6 
if a = 6 



(23) 
(24) 
(25) 

(26) 



This way of writing 7^'^^^ is rather convenient, since one can show that the 
determinants of the metric cancel, leading to a simpler final expression, de- 
pending only on O: 



Ej = - — Tr log 
2 / 27r ^ 



J+0( 



(27) 



which we use in our subsequent derivations. This equation is the A^-body 
generalization of the so called 'TGTG' formula for the Casimir interaction 



^They do contribute to the Casimir 'pressure' on each surface, though. 
^The existence of this generalization is mentioned in [7]. 
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between two bodies applied in [71 [8] (see also [9j), which in our notation 
reads: 

i?,(sW,S(^)) = ^ I ^trlog[l-G'W(^)/C(i2)(a;)G('n^)^(2i)(cc')] . (28) 

The latter is obtained from (l27j) by expanding in powers of O and summing 
up the series for the particular case N = 2: 

= - / ^ f:^tr{[0(,,)(.)(9(,,)(a;)]'=} , (29) 
k=i 

what yields (|28|) . 

Coming back to ( !27ll . we note that the form of G^""^ shall depend, im- 
plicitly, on the geometry of each surface, and in general cannot be evaluated 
exactly, except in rather simple cases. However, most properties we shall 
deal with in this section are independent of that form. 

Expression fl2711 has immediate analogues in 6? 7^ 3. Indeed, in ci = 1, we 
arrive to a result formally identical to (l27l) , after one notes that the trace only 
affects the indices that label the 'surfaces', which in this case are just points 
labelled by their coordinates x^"^ (no parameters a are involved). Besides, 
the kernel /C(u;; x) is now: 



}Ciuj-x) = (30) 

Z\UJ\ 



and G^""^ becomes: 



r -1 -1 

G^''\u)= lim . , =2|cu|, (31) 

independently of a. 

Finally, in = 2 the boundaries are curves F^") described by just one 
parameter a, and: 

^(^;x) = -^Ko(|^||x|), (32) 
Zn 

where Kq is a modified Bessel function, and G^"^ is obtained by evaluating 
the inverse of /C(aa), for which there is no general expression; we shall however 
derive its exact form for a particular case in the next section. 

Equipped with (l27l) . we can define a way to 'measure' the non-superposition 
effects. Again, we work in = 3, but the results are straightforwardly 
adapted to 6? 7^ 3: Assuming that we know , the energy corre- 

sponding to surfaces, we add an extra boundary, S(Ar+i), obtaining a new 
energy i?,({S(,)}^+^). 
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If superposition were valid, the difference between the two energies would 
be the sum of the interaction energies between and ^{a), with a = 



1, . . . , N. Thus, we introduce: 



TV 



SEj{N) ^ i?.({S(,)}tY)-i?/({S(a)}:=i)-E^K{^(^+i)'^w}) • (33) 



a=l 



Superposition is broken whenever 5Ei{N) ^ 0. Reciprocally, for the energy 
of surfaces to verify superposition we would need: 5Ei{M) = 0, for M = 



2,...,iV-l. 



The final ingredient to evaluate 5Ej{N) is obtained by applying ( l27l) to 
the + 1 surfaces. Then, we use determinant algebra to relate the determi- 
nant of the corresponding (A^ + l)*''-order matrix to an A^(*^)-order one: 



det [I^ab) + C(a6)] (^+1)^(^+1) = det [I^ab) + 



NxN 



where: 



^'(ab) = ^{ab) - 0(^a N+l)0(^N+lb) 



(34) 
(35) 



Thus, the difference between the energies for + 1 and A^ surfaces may be 
put in the form: 



f duj 




I — Tr log 




/ 27r ^ 





(36) 



where 



TV 



c=l 



{ac) 



(37) 



It is now a matter of algebra to extract the pairs interaction energy to 
show that: 

(38) 



with: 



5Ej{N) = \j ^Trlog[/ + A(u;)] 

N 

(X+0)-^](.e)0 (c N+l)0(^N+l b) 

c=l 

— 0{a N+l)0(N+l b) + C^{Ar+l a)C^(a 7V+l)^a6| 

8 



(39) 



where the / in the first factor is the identity operator on functions defined in 
parameter space (while X also acts on the indices space), the discrete indices 
are not summed, and the products are understood in the operatorial sense, 
regarding the kernels as matrix elements with continuous indices. 

In spite of the fact that the form its rather complicated, we may already 
extract some conclusions from it. The most immediate one is that for the 
strength of the non-superposition effects to be small, the magnitude of the 
matrix elements of O between the (A^ + 1)*'' surface and the previous ones 
has to be small. 

Moreover, for the correction to be smaller than the superposition terms, 
we also need (9(afe), for a, = 1, ... to be small, since these operators also 
affect the magnitude of those terms. And this is the main conclusion of this 
section, namely, that for superposition to be valid, all the matrix elements 
of O have to be small. We can see, in fact, that when that is the case, the 
form of the correction, to lowest order in the matrix elements, does depends 
on the matrix elements involving all the boundaries: 

6Ej{N) ~ -- / £ 5^tr[0(^+i,)0(,,)0(fe^+i)] . (40) 

a,b=l 

The smallness of O is what, on the other hand, renders a perturbative ex- 
pansion of the interaction energy possible. 



4 Perturbative expansion 

The condition that the matrix elements of O are small, is precisely what 
one would require in order to expand the interaction energy in powers of 
that operator. On the other hand, for O to be small, the only assumption 
available here is that the A^ surfaces S*^") are compact objects, and that the 
distance between each pair of surfaces is much bigger than the size of any 
object. Under this assumption, the norm of 0(a6) is much smaller than one, 
since the G^""^ kernel, is determined by the inverse of /C at small distances, 
while JC[ab) is, essentially, /C at long distances, and /C decreases with the 
distance. 

Excellent articles exist about the evaluation of the Casimir interaction en- 
ergy within the TGTG formula approach, by applying different expansions [TOl 
[7l[8]. We just present here an analysis of the perturbative expansion in pow- 
ers of C, from the point of view of the non-superposition effects, for the case 
of A^ boundaries. 
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The expansion yields a series for Ej: 

oo 



1=1 



where 
or: 

Ei-i = 



Em 



21 



du 
2^ 



Tr 



21 



2^ ^ 



tr 



O 



(aia2) 



(12 as) 



'lo)...0 



(atai) 



(41) 



(42) 



(43) 

This is, essentially, the long distance expansion considered in [10], although 
we only deal with the Dirichlet (strong coupling) case. 

It is worth noting at this point that the absence of explicit factors of 
the metric by no means signal a breaking of reparametrization invariance. 
Indeed, what happens is that the kernels denoted by G^""^ do have a nontrivial 
transformation properties under reparametrization, which compensate for the 
non invariance of the integrals over the parameters. 

Let us study the explicit form of the first few terms in this expansion. 
The / = 1 term vanishes, so that the lowest non-trivial order corresponds to 
/ = 2, which using ( !25i ) becomes: 



b) 



where 



a<b 

1 

~2 

X G(^)( 



duj 
2^ 



cr,(j' ,(T" ,0" 



(44) 

G^''\u;a, a') IC^'^'^u; a' , a") 
IC^'"'\u;a"',a) . (45) 



t' rr" rr"' 



To this order, the total energy is obtained as the sum of 'interaction energies' 
corresponding to the pairs, in a sort of 'superposition principle'. 
This property is violated in the next order term: 



Ei, 



du 
2^ 



E 

a,b,c 



a 



a 



X 0^ab){uj; cr, a')0(^bc){^; o-', o-")C(ca)(^; cr", cr) 



(46) 



a<b<c 
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where we have introduced a 'three-body energy interaction', E^^''^'': 

Ei-bc) ^ /^tr[G('^)(cu)/C('^')MG'('')H/C('^)(u;)G'('=)(c^)/C('='^)(a;)l . (47) 
J 2tt L J 

Incidentally, this correction coincides with ( 1401 ) when one considers + 1 
surfaces, as it should be, since on should expect that the lowest order violation 
to the non-superposition comes from the lowest non-quadratic term in the 
energy. 

A fundamental ingredient in the calculation of the different terms in the 
expansion for Ej is the kernel G^"^(a;; a, a'). The form of that kernel depends 
strongly on the number of spatial dimensions as well as on the shape of the 
surface itself. Universal statements can only be made if more assumptions 
about the surfaces are made. However, based on the same assumption used 
to perform the series expansion, we may simplify the previous expressions 
further. Indeed, denoting by x^^^ the barycenter of the T.^"''^ surface, we can, 
in the expressions above, use the approximation: 

IC'^''''\u;a',a") ~ ^(a;; x^^^ - x^^^) . (48) 

This is justified by the following reason: we are assuming that |x('^'-x(^)| >> 
Ri^)^Rib)^ where R^"^ denotes the minimum radius for a sphere 5**^"^ cen- 
tered at x^"-*, which encloses E^. Then we may replace y^"''(o") x.^'^^ and 
yW(cr') — > x^''^ since K.^'^^^u; a' , a") is (under the previous assumptions) 
approximately constant inside S^^\ 

Using this approximation inside the expression for E^°-'^\ we see that it 
may be written as follows: 

E^'^'^ ^ Jd'xJ d'yp^'^\u;x)Viu;x-y)p^''\y) , (49) 

where we introduced: 

Qaiuj) = J d'^a J d^a' G^''\io;a,a') (50) 

and 

V(uj; X - y) = -i [k{uj; x - y)] ' . (51) 

Thus, at this order, we see that the interaction energy for the a, b pair, may 
be regarded as arising from integral over u of the interaction energy for a set 
of pointlike charges located at x^^^ and x'^''^, whose strengths qa{i^) and qb{i^) 
are determined by the geometry of the respective surface. 
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On the other hand, the explicit form of the interaction potential is: 

1 p-'iH\^-y\ 

Vi.;.-y) ^ (52) 

hence, the interaction is always attractive. The integrals over x and y have 
been used in order to make it clear that each surface behaves as a sort of 
point-like charge. Of course, the same approximation may be used to simplify 
the form of the higher order terms. 

The form of G^""^ is not known exactly in general, except for particular 
situations, like the d = 1 case, which we consider now: 



4.1 d=l 

As a first test, we consider the case of two mirrors in 1 + 1 dimensions. The 
operator O is just an cu-depending matrix, with matrix elements 0(^ab){^) = 
g-|i^||a;('')-a;('')|. exponential decay assures the convergence of the perturba- 
tive expansion, regardless of the relative distances between the mirrors. 

In this situation, the first (superposition) expression for the energy cor- 
responding to two point-like objects (mirrors) located at x^^^ and x*^^^ adopts 
the form: 

= j ^(2|u;|)V(c.^;x«-x(2)) , (53) 

where 

n.;x.«-x.^)) = --^^^ (54) 
Assuming that the distance between the mirrors is a, we see that. 

^(12^ 1 /•+°° du ouAn 1 0.07958 , , 

to be compared with the exact result: ^ = ~2fa — which is bigger 

by approximately a sixty percent. 

It is possible to calculate, for this case, all the higher order corrections 
exactly; only the even orders yield non-vanishing contributions, which are 
given by: 

""''^^ = -21 y_ 2^^' ' = -WiTa- ^'^^ 
Then one sees that their sum: 

El 1 In TT 

1=1 1=1 
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which is the exact result. 

Besides, when more than two mirrors are considered, the energy becomes 
equal to the sum of the Casimir energies corresponding to the pairs formed 
by neighboring mirrors: 



N~l 



TT X ^ 1 



a=l 



4.2 d = 2 

For a radius R circle, using the angle as parameter, we find: 

For an infinite line, parametrized by a G {—oo, +oo), the result is instead: 

G^''\u;a,a') = 2 { - ^ + u;')Ko{Ma - a'\) . (60) 

It is straightforward to check that, in both cases, one is in a situation of an 
O with small norm. Thus, in the perturbative expansion, when the surfaces 
are bounded and very far away, using the approximation that follows from: 

JC^ab){uj; Ix^'^) - x^")]) ~ IC{u; jx^'^) - x(^)|) (61) 

we may obtain an approximate expression for the case of circles. Denoting 
by r]^""^ and x*^"^) the center and radius of each circle, and assuming that 
|x('^) — x*^*)] >> max{r]^"-^}, we have for the pair interaction energy: 



2 J 2n /o(|^|?7('^))Ko(kl?7('^))/o(kl^('))^o(k|r/W) 

For the case of just two circles, 1 and 2, say, one can sum the series corre- 
sponding to the different powers of O. The result is: 



~ -2 / ^ 



1 - 



[Ko{\co 


|xW - X 


W|)]2 




/o(|^ 


?7("))A'o(|u;|? 


7W)/o(|a; 


r/W)i^o(kl??(')) 



(63) 
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4.3 d = 3 



We now deal with the case of surfaces in d = 3. The long distance approx- 
imation requires the evaluation of the integral of G^""^ over the parameters; 
for a case of a sphere, that integral is: 

This is smaller than ]C(ab) for distant surfaces, as a straightforward test shows. 

If one assumes instead that the surfaces are really composed of small, 
weakly coupled surface elements, we may in fact use local approximations 
for the G^'^^ kernels. In this local approximation means that the 

kernel is concentrated around a = a': 

G^"^) {uj; a, a') ~ r/^"^) {u, a) 5^^^ {a - a') . (65) 

where a = 1,2, and rj^'^'^ will be determined now: Recalling that G^""^ is 
defined as the inverse of /C(aa)(a;; a, a'), we explore the form the latter in the 
neighborhood of a given point in the surface S^, the one characterized by 
the parameter a: y^^\a). Close to that point, we derive the approximate 
expression: 



£^ gik||-0<,y(")(<7)K-<7'-) 



1 



exp[-\u\Jgl%)6a-6(T(^] ,(66) 



where ky is the projection of the momentum along the tangent plane at the 
point y(")(f^), and Sa"" = a" - a'". 

In the assumption that there is no appreciable momentum flux between 
the different surface elements, we end up with the expression: 

This yields: 

G^''\uj;a,a') ~ 2^ g(-){a)\u\ 5^^\a - a') ^ r]{u, a) = 2\uj\ g^^^a) . (68) 

Let us first assume that we have two surfaces, S*^^) and T,^^\ and con- 
sider the second order expression for the interaction energy, using the local 
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approximation for the kernels G^^^ and G'^^^ We see that their interaction 
energy at this order becomes: 

(69) 

which, performing the integration over u, results in the following expression: 
= J d^a^gWia) J d^a' ^ g(^){a')V{a, a') (70) 
where: ^ ^ 

^^""'"^^ " "32^ |yW(^^) -y^'H^OP ^'^"^^ 

which looks like a kind of local-potential, van der Waals like interaction. 

Finally, let us consider the case of infinite parallel plates, within the 
quadratic approximation, using two different approaches. Obviously, in this 
case, the planes cannot be regarded as small surfaces and, even though the 
superposition approximation may be valid, certainly the planes cannot be 
regarded as point-like objects. 

It is clear that, using as parameters the coordinates Xy"^ on each mirror, 
we have: 

G^^\u- xj'^), yff)) = / H 2 e-ii-(-r-r) , (72) 

while for O the result is: 

0,.,(.;x<"',y!") = / |i ^-V^F^M->-.«,«.,.(«i">-.i-') _ ,73) 

where z^"-^ is the position (on the third axis) of each plane. We see that, even 
in this case, the norm of the operator is small. 

Indeed, inserting this into (l45ll , we get for Sq, the energy per unit area: 

1 f dU f d^k\\ -2a^m^ 1 ^r.^oo -3 /^.N 

^» = -2 / ./ (2^? ^ ^" =i6A5"-°-°°«="' ■ f^"' 

to be compared with the exact result, that is £i-2 — —0.006854 a"'^, what is 
a signal that the corrections are small. 

If, on the other hand, we imagine each mirror as composed of weakly 
interacting infinitesimal surface elements (not a conductor), and apply the 
superposition result to a system composed of all the surface elements, then 
the energy per unit area to the first non-trivial order, £1-2, may be obtained 
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by integrating the interaction energy between a single point on a mirror and 
all the points in the other. This corresponds to the following integral: 




The integral can be evaluated exactly, yielding: 

8,,, = ^ -0.00422 a-^ (76) 

which is different than the previously obtained result, as it corresponds to a 
different material. 

5 Conclusions 

We have obtained an expression that measures the departure from super- 
position in the interaction Casimir energy corresponding due to a massless 
scalar field in the presence of > 2 Dirichlet surfaces. We have found that 
the most general condition under which the non superposition effects can be 
regarded as small corresponds to a number of small surfaces separated by 
long distances. Under this assumption, one may construct a perturbative 
expansion, as a series in the operator O. 

The condition on that operator manifests itself in a different fashion, 
depending on the number of spatial dimensions, hid = 1, since the size of the 
mirrors is zero, one is in the best possible situation, namely, the perturbative 
expansion is always reliable. 

In d = 2 and d = 3, on the other hand, one can always obtain con- 
ditions under which the expansion should be reliable (although the rate of 
convergence depends on d). 

An interesting conclusion one can extract is that the knowledge of the 
interaction energy for surfaces is useful to calculate the one for A^ + 1 
surfaces only when all the surfaces are widely separated. 
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